Finite-difference frequency-domain

The finite-difference frequency-domain (FDFD) is a numerical solution for problems usually in electromagnetism and sometimes in acoustics, based on finite-difference approximations of the derivative operators in the differential equation being solved.

While "FDFD" is a generic term describing all frequency-domain finite-difference methods, the title seems to mostly describe the method as applied to scattering problems. The method shares many similarities to the finite-difference time-domain (FDTD) method so much of the literature on FDTD can be directly applied. The method works by transforming Maxwell's equations (or other partial differential equation) for sources and fields at a constant frequency into matrix form Ax = b. The matrix A is derived from the wave equation operator, the column vector x contains the field components, and the column vector b describes the source. The method is capable of incorporating anisotropic materials, but off-diagonal components of the tensor require special treatment.

Strictly speaking, there are at least two categories of "frequency-domain" problems in electromagnetism.[1] One is to find the response to a current density J with a constant frequency ω, i.e. of the form \mathbf{J}(\mathbf{x}) e^{i\omega t}, or a similar time-harmonic source. This frequency-domain response problem leads to an Ax = b system of linear equations as described above. An early description of a frequency-domain response FDTD method to solve scattering problems was published by Christ and Hartnagel (1987)[2]. Another is to find the normal modes of a structure (e.g. a waveguide) in the absence of sources: in this case the frequency ω is itself a variable, and one obtains an eigenproblem Ax = \lambda x (usually, the eigenvalue λ is ω2). An early description of an FDTD method to solve electromagnetic eigenproblems was published by Albani and Bernardi (1974).[3]

Contents

Tips for Implementing the Method

  1. Use a Yee grid because it offers the following benefits: (1) it implicitly satisfies the zero divergence conditions to avoid spurious solutions, (2) it naturally handles physical boundary conditions, and (3) it provides a very elegant and compact way of approximating the curl equations with finite-differences.
  2. Much of the literature on finite-difference time-domain (FDTD) applies to FDFD, particularly topics on how to represent materials and devices on a Yee grid.

Literature

References

  1. ^ J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade (2008). Princeton Univ. Press. ed. Photonic Crystals: Molding the Flow of Light, 2nd edition. pp. 688–696. http://ab-initio.mit.edu/book/. 
  2. ^ Andreas Christ and Hans L. Hartnagel (1987). "Three-Dimensional Finite-Difference Method for the Analysis of Microwave-Device Embedding". IEEE Transactions on Microwave Theory and Techniques 35 (8): 688–696. Bibcode 1987ITMTT..35..688C. doi:10.1109/TMTT.1987.1133733. 
  3. ^ M. Albani and P. Bernardi (1974). "A numerical method based on the discretization of Maxwell equations in integral form". IEEE Transactions on Microwave Theory and Techniques 22 (4): 446–450. Bibcode 1974ITMTT..22..446A. doi:10.1109/TMTT.1974.1128246. http://ieeexplore.ieee.org/xpl/freeabs_all.jsp?arnumber=1128246. 

See also

Finite-difference time-domain method